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An introduction into the theory of boundary critical phenomena and the applica- 
tion of the field-theoretical renormalization group method to these is given. The 
emphasis is on a discussion of surface critical behavior at bulk critical points of 
magnets, binary alloys, and fluids. Yet a multitude of related phenomena are men- 
tioned. The most important distinct surface universality classes that may occur 
for a given universality class of bulk critical behavior are described, and the re- 
spective boundary conditions of the associated field theories are discussed. The 
short-distance singularities of the order-parameter profile in the diverse asymptotic 
regimes are surveyed. 



1 Introduction 



Until the end of the 70ies at least the subject of boundary critical phenomen; 




attracted only rather limited attention, even though a number of pioneering 
papers^ had already been written on it. This was due to several reasons. 
Perhaps the most important one was the complete lack of sufficiently precise 
experimental work. Such work even seemed beyond reach of the then available 
experimental possibilities. Another reason was that the powerful machinery 
of the field-theoretical renormalization group (RG) approach had not yet been 
extended to systems with boundaries, and hence could be utilized for system- 
atic studies of boundary critical phenomenaTli3 only later. Finally, also the 
seminal work of Belavin, Polyakov, and Zamolodchiko^a on two-dimensional 
conformal field theories had yet to be performed and extended to systems with 



Today the situation has changed considerably. The subject has become 
a very active research field, attracting many scientists with a broad range of 
diverse backgrounds. This includes the experimental and theoretical condensed 

tThis article is based on an invited talk presented at the Third International Conference 
"Renormalization Group — 96", Dubna, Russia, August 26 — September 1, 1996. 
"Due to lack of space and the great number of papers concerned, it will not be possible to 
present here an adequate account of these earlier contributions and their importance. For 
the same reasons I shall not in general be able to refer to many original papers. My citattmi 
philosophy will be to refer, wherever possible, to appropriate review articles (such asBEl) 
by whose extensive lists of references the interested reader can easily trace the literature. I 
apologize to all colleagues whose papers could not explicitly be cited here. 




matter physicist wha_is interested in the critical behavior of magnets and alloys 
with free surfacefflolij just as well as the physical chemist investigating the 
critical adsorption of binary fluids on walls and interfacesjLj or the theoretical 
high-energy physicist studying field theories on manifolds with boundariesEZI 
Further, there is a wealth of related phenomena. The reader can get a first 
impression from the following incomplete list of topics belonging to the realm 
of boundary critical phenomena: 

m 
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(i) critical behavior of magnets and alloys with free surfaces}! 

(ii) critical adsorption of fluids on walls and interfaces ^0 

(iii) critical adsorption of polymers on wallsfH 

(iv) percolation transitions in systems with boundaries^! 

(v) quantum impurity problems (such as the Kondo effect)[' 

(vi) Schrodinger representation in renormalizable field theories^ 

(vii) Casimir effect@ 

(viii) relaxation processes starting from initial non-equilibrium states0@ 
(ix) dynamical critical phenomena in systems with boundaries l'E3 El 

The common feature of these phenomena is that their physics on large 
length and time scales is described by nearly critical (or 'massless') field theo- 
ries on manifolds with boundaries. These boundaries may be genuine surfaces 
[as in examples (i), (iv), (vii), (ix)], walls [as in (ii), (iii), (vii)], or interfaces 
separating a nearly critical phase (such as a binary mixed fluid near its conso- 
lute point) from a non-critical spectator phase (such as vapor), where the latter 
example applies to the case_of critical adsorption (ii) of a biliary fluid mixture 
at its critical end pointtijE3 In quantum impurity problems [topic (v)] , the 
role of the boundaries is played by the impurities, which turn out to be equiv- 
alent to boundary conditions of the resulting effective 1+1 dimensional field 
theories. In the examples (vi) and (viii) one is dealing with the time evolution 
from a given initial state at time f^; here the space-time hyperplane t — t{ is 
the analog of a boundaryEjo 

The purpose of this paper is to give a brief introduction into the field of 
boundary critical phenomena and to survey some recent pertinent RG results. 
The emphasis will be laid on a discussion of topics (i) and (ii). This also serves 
to explain the basic questions arising in, as well as general aspects of, boundary 
critical phenomena. In addition, typical theoretical strategies will be described 
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and some representative experimental results mentioned. A detailed exposition 
of the other topics, (iii)-(ix), is beyond the scope of the present article. The 
reader should consult the cited literature for more information about these. 

2 Surface effects 

Owing to their conceptual simplicity, ferromagnets are a convenient starting 
point. In studies of bulk critical behavior, two important properties of real fer- 
romagnets are usually ignored: their finite size and the presence of surfaces. 
To gain insight into the significance of these properties, consider a ferromagnet 
of finite linear extension L, where L is large in comparison with the charac- 
teristic microscopic scale, the lattice constant a. For concreteness, we take 
a d dimensional cube. Its 2d faces make up the boundary of total surface 
area A — 2dL d ~ 1 . We presume that the microscopic interactions are of short 
range. As a microscopic model one can imagine an Ising model on the lattice 
[0, L] d n Z d , with a Boltzmann factor exp(— Tii a t) of the form 

Wlat = \ Kij Si Sj . (1) 

Here i labels the sites, and sa — ±1 are spin variables. Away from the surface, 
the interaction constants are translation invariant, K^j = K(i—j), but close to 
it K^j will in general also depend on the distance of the bond's center from the 
surface. Restricting ourselves to nearest-neighbor (nn) interactions, we assume 
that K^j takes the value K\ for all nn bonds between two surface sites, the 
'bulk value' K for all other nn bonds, and is zero otherwise. It is understood 
that both K and K\ are ferromagnetic (< 0). No restriction is imposed on the 
boundary variables (free boundary conditions ).[] 

The presence of surfaces entails that thermal averages of local densities 
deviate from their values deep inside the sample. Far away from the critical 
point K — K c , this disturbance can penetrate into the sample only up to a 
distance comparable with the interaction range, for the correlation length £ is 
of the same (microscopic) order. However, as the temperature T approaches 
the bulk critical temperature T c , £ grows, getting macroscopically large. Hence 
the boundary region affected by the surface acquires a thickness of the same 
macroscopic order. An obvious first consequence is: 

b If required, pure bulk phases with up or down magnetization are selected by means of a 
homogeneous magnetic field H which approaches ±0 after the thermodynamic limit has been 
taken. 
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• Local densities such as the order parameter m(x) or the energy density 
e(x) become inhomogeneous on the scale of £. 

In our case, m(x) and e(x) correspond to (coarse-grained versions of) the local 
magnetization m{i) = (sa) and energy 

£ (i) = E ( s i s j)> ( 2 ) 

mi bonds of t 

respectively. 

The position-dependence of the magnetization for T <T C along a line per- 
pendicular to one pair of faces, and sufficiently far away from the others, is 
shown schematically in Fig. 1 for the case £ <C L. Deep inside the bulk, m(x) 




Figure 1: Schematic order-parameter profile m(z) between two parallel faces, where z is the 
distance from the left face. Up to exponentially small correction terms due to the faces not 
shown, the profile is independent of the coordinates perpendicular to z. 

agrees with the bulk value m& up to corrections ~ exp(— z/£) (and analogous 
ones with z replaced by the distance to the other faces). Owing to the re- 
duced coordination number at the surface, one expects the profile to decrease 
upon approaching the surface provided the interactions are not too strongly 
enhanced near the surface. However, it is also possible that the magnetization 
profile bends upwards near the surface. This may happen when the enhance- 
ment of surface interactions is sufficiently strong, or when the magnetization is 
locally increased by a magnetic field acting only on spins in the vicinity of the 
surface. As we shall see later, there may even be changeovers between different 
asymptotic forms of large-scale behavior in the region z < £ if certain surface- 
related crossover lengths become much larger than the lattice constant a. Of 
course, on the scale of a, the behavior of the densities depends on microscopic 
details (chosen interactions, lattice type, etc). 
A second, equally obvious, consequence is 
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• the appearance of surface corrections in integrated densities. 



For example, the total order parameter per volume V (= L d ) can be written 

as 

If A 

— / d d x m(x) w m b + — m s + . . . (3) 
v Jm v 

in the limit L — > oo. Here M. denotes the region C M d occupied by the system. 
The ellipsis stands for (edge) contributions ~ L~ 2 and faster decaying ones. 
The quantity m s is called surface excess magnetization. An analogous decom- 
position into a bulk term, and an excess contribution / s , called surface free 
energy density, holds for the reduced free energy 

= -lnTre~ Wlat (4) 



k B T 

and other quantities. 

Surface corrections like AV~ 1 f s are small for large L as long as £ 
For temperatures T sufficiently close to the pseudo-critical temperatur 
Tl = T c + 0(L~ U ) of the finite system, the correlation length is of the same 
order as L. Then the separation (||) into bulk and surface contributions looses 
its meaning. Finite-size effects become important. These have been surveyed 
by Dohm recentljO. Since our subject here is boundary effects, we will always 
assume that (ci. Hence we first take the thermodynamic limit L — > oo and 
only then consider the approach to the critical point. By this procedure finite- 
size effects are eliminated. It amounts to the study of semi- infinite systems. 
In most of the following we will restrict our attention to these. 

The critical behavior of bulk densities is described by familiar power laws. 
For the order parameter mb, and the singular parts of the energy density Eb 
and free energy fb, one has, in zero magnetic field, 

m 6 wM_|r| /3 asr^CT, (5) 

and 

ef^^M 1 -, Jf ng) « B ± |t| 2 - asr^Oi, (6) 

respectively. Here a and /3 are two independent standard critical exponents in 
terms of which all other commonly introduced bulk critical exponents can be 
expressed]^ Af_ and B± are nonuniversal (metric) factors. 

c We assume that the bulk dimension d is between the upper and lower critical dimensions 
d* and d t , so that hyperscaling is valid. Thus a = 2 — du and j3 = (u/2)(d — 2 + r]), where 
v and r) are the usual correlation-length and correlation exponents, respectively. 
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It is natural to expect that local densities, taken at points within a distance 
z <C £ from the surface, also display critical behavior characterized by power- 
law singularities, as the bulk critical point is approached. In other words, a 
third consequence should be 



the appearance of new critical indices characterizing the critical behavior 
of surface quantities. 



Thus, for the local surface order parameter mi = m(z = 0), one anticipates 
that 

mi ~ \t\ Pi as t -> 0" , (7) 

where ft in general should be different from ft 

This is indeed the case. It sufficesJaere to mention just a few exemplary 
sources of evidence. From exact worlc3 on the semi-infinite Ising model in 
al = 2 dimensions it is known that ft = 1/2, whiphJs. to be compared with 
the Onsager value (3 — 1/8. Even mean- field theoryMElfl gives different values, 
namely ft = 1 and (3 = 1/2, independent of d. The first clear experimental 
evidence for a value of ft ^ (3, to my knowledge, is due to Alvarado et ale3 
These authors investigated the temperature dependence of mi at the (100) 
surface of a Ni ferromagnet via spin-polarized low-energy electron diffraction. 
They found ft = 0.825^;^ , a value much larger than the bulk exponent (3 rs 
0.36 of the three-dimensional 0(3) Heisenberg model, with which the observed 
bulk critical behavior agrees quite well— Their result for ft is in reasonable 
agreement with the theoretical estimated ft (d = 3, n — 3) = 0.84 ± 0.01, which 
was abtained by combining the results of. an e = 4 — d expansion to second 
ordeflo with those of a d — 2 expansion^ to the same order. (Here n stands 
for the number of components of the order parameter.) 

More recently, the method of X-ray scattering under grazing incidence has 
been employed with impressive success to study the surface critical behavior 
of the binary_alk|y FeAlallj As suggested in the theoretical work of Dietrich 
and Wagnerpj c3 this technique lends itself well to accurate investigations of 
surface critical exponents. It enables one to determine several surface critical 
exponents independently, so that scaling relations can be checked. 

EeAl has a rather rich phase diagramtj The measurements of Mailander 
et aoii3 were performed at the continuous bulk phase transition between the 
high-temperature phase with B2 ordering and the low-temperature phase with 
DO3 ordering. This transition involves a two-component order parameterpl 
its bulk exponent (3 should agree with the value (3(3, 2) « 0.345 of the 0(2) 
|(/)| 4 model in d = 3 dimensions.^] This appears to be the case, although the 

d The Hamiltonian of an appropriate continuum field theory for this bulk transition clearly 
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precision of thejsxperimental estimates is not sufficient to discriminate between 
the Ising valued /3i(3, 1) w 0.33 and the one for n — 2_given above. In any case, 
a much larger value of pi was found again, namelja pi = 0.75 ± 0.02. This 
result is close to the estimatel ft. (3, 2) w 0.8 obtained from the e expansiontKij 

Monte Carlo calculations for the Ising cas£lrEa also gave clear evidence 
for values of Pi ^ p. 

Summing up, it may be said that there is indisputable evidence for the 
fact that the surface exponent pi generally is different from the bulk exponent 
p. However, this is only part of the story. A central element of the modern 
theory of bulk critical phenomena is the division into (bulk) universality classes. 
An obvious question to be asked is whether a similar classification can be 
accomplished for surface critical behavior at bulk critical points. The answer 
is yes. We call such classes surface universality classes. As it turns out, 

• for a given bulk universality class, there exist in general several distinct 
surface universality classes. 

To explain this we must consider specific models. Instead of continuing to work 
with lattice models such as the one introduced in Eq. (Q), we will now turn 
directly to appropriate continuum models describing the large-scale physics. 

3 Surface universality classes 

We presume that on large scales no long-range interactions must be taken into 
account. This requires that the following necessary conditions are fulfilled: 

(I) All pair interactions must be of short range or at most have long-range 
parts that are irrelevant in the RG sense (anywhere, i.e., both in the bulk 
and in the vicinity of the surface) . 

(II) All boundary-induced contributions to the interactions must be of short 
range in the distance z from the surface or at most have parts of long 
range in z that are irrelevant in the RG sense. This must hold, in par- 
ticular, for any one-body potential associated with the boundary. 

One-body potentials decaying as a power as z — > oo occur in the case of 
fluids bounded by walls. However, the exponents ui usually are larger than the 
critical value u c = P/v below which such long-range tails must be expected to 
become relevant (see, e.g., pages 210-213 of Ref. §). 

will not be totally 0(2) symmetric. Wawever, symmetry-breaking terms such as cubic 
anisotropics are believed to be irrelevant^ for n = 2. 
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The large-scale physics of systems satisfying (I) and (II) (which the lattice 
model of the previous section evidently does) may be expected to be described 
by a continuum field theory with a local action of the form 

U= I dVC{x)+ I dA£i(x) , (8) 

where C(x) and C\(x) are functions of the order-parameter field 4>(x) and its 
derivatives. For a general d dimensional manifold A4 with boundary QA4 = B, 
the representation-invariant volume and area elements are given bjtj dV ~ 
^/gd d x and dA = ^fgs d d ~ 1 x\\ , respectively, where fl'(x) is the determinant 
of the metric tensor while gs is its analog for the induced metric on B. Wc 
will restrict ourselves to manifolds M. C M d in the sequel, and unless stated 
otherwise, we will also not consider curved boundaries. Working with semi- 
infinite systems, our standard choice of M. will be the d dimensional half-space 
R d + = {(x\\,z) e M d | < z < oo}, withS given by the z = plane. Thus dV 
and dA can simply be read as d d x and d d ~ 1 x\\, respectively. 

The 'bulk density' C(x) must be chosen in such a way that a proper de- 
scription of bulk critical behavior results. In the case of a usual critical point 
of systems belonging to the universality class of the isotropic n-vector model, 
this tells us to choose 

C{x) = \[\/cp{xf +W[0(x)] (9) 

with 

U{4>) = |r o 2 + ^oH 4 , (10) 

up to symmetry-breaking terms. 

In the overwhelming part of the following it will not be necessary to include 
derivative terms in the 'boundary density' L\. Hence we choose it to be of the 
form 

£i(x) =Ui[(j>{x)] . (11) 

In fact, for the Hamiltonians representing the surface universality classes we 
will be concerned with, the choice 

Ui(4>) = ic o 2 -/ii, o (12) 

is general enough. To understand why no |0| 4 boundary term is included, 
one should note that the associated coupling constant would have momentum 
dimension e— 1. Thus at least for sufficiently small e = 4— d, such a term should 
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be irrelevant^ A boundary contribution of the form (f>d n 4> cannot be excluded 
on the basis of power counting alone since it involves a dimensionless coupling 
constant. However, it can be ruled out on other grounds (it is redundant )ljErLil 
The phase diagram of the semi-infinite n- vector model defined by Eqs. (|8|)- 
( |l2| ) with hi Q = is shown in Fig. 2. The vertical axis corresponds to Tq; 




surface enhancement — c — > 
Figure 2: Phase diagram of the semi-infinite n-vector model for d — 1 > d*(n). 



in identifying it with T, we have made use of the fact that tq is a linear 
function of T near T c . The horizontal axis corresponds to cq, whose negative 
is a measure of the enhancement of the surface interactions. The variable c is 
proportional to cq — c sp , where c sp is chosen such that the point Sp is located at 
c = 0. The lines labeled "ordinary", "extraordinary", and "surface!' represent 
continuous phase transitions that have been given these namestlj The lines 
meet at a multicritical point Sp, which describes the "special" transition. The 
ordinary transition occurs for subcritical surface enhancement c > 0; it is a 
transition from a surface-disordered (SD), bulk-disordered (BD) phase to a 
surface-ordered (SO), bulk-disordered phase. On lowering the temperature at 
fixed supercritical enhancement c < 0, one first enters a surface-ordered, but 
bulk-disordered, phase at T — T s (c). This so-called "surface" transition is like 
a bulk transition of a d — 1 dimensional system. (The bulk correlation length 
of the d dimensional semi-infinite system remains finite on the transition line, 
as long as c stays away from its critical value zero.) As T is lowered further at 



e A |</>| 4 boundarjj— tmn is needed to describe certain surface universality classes of bulk 
tricritical system Mi 11 ! in dimensions d < 3. Of course, in this case a bulk term is 
required as well. 
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fixed c < 0, the extraordinary transition takes place. 

At bulk criticality, T = T c , we thus have three distinct types of surface 
transitions: the ordinary, the special, and the extraordinary one, provided d 
exceeds their respective lower critical dimensions (led). The led of the surface, 
extraordinary, and special transitions is 2 for n = 1, but 3 for n > l.This is be- 
cause a SO/BD phase cannot appear (except for infinite surface enhancement), 
unless the surface dimension d— 1 is larger than d*(n) = 2 — 5 n i. Hence, for the 
isotropic Heisenberg case n = 3, only the ordinary transition remains in three 
dimensions .£] In the particular case d— 1 = n = 2, a phase with quasi- long- 
range order at the surface is possible. This possibility will not be considered 
here. Note also that the limit n — ► is known to describe the adsorption of 
polymers on a wallE3 For this problem an analog of the multicritical point Sp, 
corresponding to an adsorption threshold, exists for d > 2. 

For given d and n, each one of the above-mentioned three types of transi- 
tions of our model is characteristic of a separate surface universality class, also 
called ordinary, special, and extraordinary, and described by (h^o = 0) RG 
fixed points 

^ord : C ord = 00 > Kp : C sp = I Kx : C *x = ~°° > ( 13 ) 

respectively, with the indicated fixed-point values of c. 

Let us briefly summarize some important properties of these transitions. 

3.1 The ordinary transition 

At the ordinary transition there is a single relevant surface scaling field gi, 
which varies ~ h\ o for small hi q. All conventionally defined surface critical 
exponents can be expressed in terms of its scaling index A° rd /i/ and two in- 
dependent bulk exponents. For example, the critical surface pair correlation 
function, the singular part of the surface energy density, and mi behave as 

(d>(x h 0) ( f>(O,0)) T=Tc ^\x ]l \-^- 2+ ^ , (14) 

4 8ing) «S ± |r| 2 - Q , (15) 

and 

mi « \ t\^ a±( ffl |r|- A ? rd ) , (16) 

■^A SO/BD phase becomes possible in the Heisenberg case for d > 2 when there is an easy- 
axis spin anisotropy at the surface with a sufficiently enhanced surface interaction constant. 
For a critical value of the enhancement the transition temperature of the Ising-like surface 
transition coipjides with the T c of the Heisenberg bulk transition. This defines an anisotropic 
special point£j which must not be confused with the isotropic ones, Sp, we have in mind 
here. 
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with 

/9i = ^(d-2 + »jy) = (d-l)i/-Ai. (17) 

The exponent values of /3i given in the previous section all referred to t.hr. 
ordinary transition. The agreement of the field-theory estimatefl'ErEHi3'C§E3 
with those obtained by other means is quite—good, though one is far away 
from the precision achieved in the bulk casco because the field-theoretical 
calculations have been carried out only to two-loop order so far. For the Ising 
case, one finds /3° rd (3, 1) « 0.8. Comparisons with recent Monte Carlo and 



other estimates may be found in Refs. 51. 52, and 56. Experiments with 
binary liquid mixtures yielded the estimattEl /3° rd = 0.83 ± 0.05. (For another 
recent experiment, see Rcf. |58| .) 

Sincc_scaling laws such as ( J17j ) can be derived quite generally from the RG 
equations™, they must hold to any order of the e expansion (and the d — 2 
expansion^, for n > 2). From the RG equations one can also conclude that 
scaling functions such as a± are universal up to a fixing of scales. In the bulk 
case there are two nonuniversal metric factors to be fixed ( "two-scale- factor 
universality"). Here there is one additional nonuniversal factor, associated 
with the relevant surface scaling field g\. Hence we have a (2 + 1)- scale- factor 
universalityu 

To explain another important property in an elementary fashion, we note 
that the order parameter profile takes the scaling form 

m(z,r)^M^\rfP(z/0 with P(C) « { J^Cft-ffl/* ™ £ 2 ( 18 ) 

where the short-distance singularity of -P(C) follows from the requirement that 
m(z <C £) has the temperature dependence r^ 1 of mi. Since /3f rd > f3, P{0 + ) 
vanishes. That is, the order parameter satisfies an asymptotic Dirichlet bound- 
ary condition. This argument can be generalized in a straightforward fashion 
to multi-point correlation functions, both for r > as well as for r < 0. As 
first shown in Ref. ^, the short-distance singularity <f> z ^o ~ z^-M v can be 
systematically obtained by means of a short-distance expansion in terms of 
boundary operators (cf. Refs. ^ and [Tt], and Sect. 4 below). Explicit results 
for the scaling function P(Q obtained by RG-improved perturbation theory to 
one-loop order may be found in Ref. ^|. 

3.2 The special transition 

This transition involves besides the analog of gi, which we call h±, a second 
relevant surface scaling field, c. One is led to scaling forms such as 

mi(r,c) w |T|^±(cM-* > fci|T|- A i P ) ■ (19) 
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Standard arguments based on crossover scaling imply that the crossover expo- 
nent $ governs the behavior of the line of surface transitions near Sp, giving 
[T s (c)~T c ]/T c ~ |c|V*. 

In mean-field theory /3^ p = $ = 1/2, independent of n. The e expansion is 
knoven to second order. Setting e = 1 in the corresponding 0(e 2 ) expressions 
gives /?i P (3, J4 i 2 ^ 0.22... 0.25, in reasonable agreement with various Monte 
Carlo resultspM3 which yielded values in the range 0.19 . . .0.24 (even though 
the error bars reported in some of these works arejiiuch smaller). The agree- 
ment is less satisfactory for ^—Earlier simulations!!! gave $(3, 1) w 0.59, but 
more recent Monte Carlo worlc3 suggests a value < 0.5. The discrepancy with 
the reported e expansion estimate is presumably due to the unusually large 
0(e 2 ) term of $. This view is supported by the results of a recent paperfiS in 
which the massive RG approach for fixed dimension d = 3 has been extended 
to semi-infinite systems. From a two-loop calculation and subsequent Pade- 
Borel analysis the value $(3, 1) w 0.54 was found. If such elaborate techniques 
are used to extrapolate the e expansion to d = 3, one can get down to similarly 
small estimates of $. To significantly improve the accuracy of such field-theory 
estimates, calculations to higher loop orders are required. 

The scaling form ([l8]) of the order-parameter profile carries over to the 
present c = case (with n — 1), but now the exponent /3l P — f3 of the short- 
distance singularity is negative for 3 < d < 4. Hence the analog of -P(C) tends 
to +oo as £ — ► 0. Only when /3^ p = (3 (as in mean-field theory and hence for 
d > 4), do we have a Neumann boundary condition d n (f> = 0. Again, these 
scaling considerations can be extended to TV-point correlation functions, and 
the short-distance singularities can be derived in a systematic manner via a 
boundary-operator expansion (cf. Sect. 4). 

3.3 Extraordinary and normal transitions and critical adsorption 

The special feature of the extraordinary transition is that there is spontaneous 
symmetry breaking and long-range surface order both below and above T c . Bray 
and Moorio argued that the important point is just that mi remains nonzero 
at the transition. Instead of having long-range surface order caused by spon- 
taneous symmetry breaking for T > T c , one could as well have surface order 
induced by a surface magnetic field hi^. In other words, the surface critical 
behavior of the extraordinary transition (supercritical enhancement c < and 
^1,0 = 0) should be the same as in the case of subcritical enhancement c > 
with hifi > 0. Applying scaling arguments to the latter case, they predicted 
the behavior m^ smg ' ) ~ |r| 2 ~", corresponding to /3J X = 2 — a. 

The case c > 0, hi t o ^ is quite normal for bounded fluids and binary 
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fluid mixtures in contact with a wall, since for these a surface ordering field 
(and other terms breaking the Z12 symmetry) generically should be present 
even at bulk coexistence. For this reason, the transition with 7^ and 
arbitrary c < 00 has been termed "normal" Ocj Thus the claim is that the 
extraordinary and normal transition are representative of one and the same 
surface universality class, provided both are possible. (Note that the Led. of 
the normal transition is d — 2, whereas it is d = 3 for the n = 1 extraordinary 
transition.)^ 

In subsequent work by Burkhardt and the present authoS it has been 
possible to demonstrate the extraordinary-normal equivalence for the Ising 
case in an exact manner. Likewise, the asymptotic behavior at the normal and 
extraordinary transitions, 

mi -m ( ™ K) «S±|t| 2 - q (20) 

with 

m^ rcg) = mi, c + A\t + A 2 t 2 + ... , (21) 

is well establishedBEH The ratios of the coefficients B±, B± and B± of 
the |t| 2-q singularities in Eqs. (|J), (jl5| ), and ( |20| ) are all given by the same 
universal bulk ratio, i.e., B + /B^ = B + /B- = B+/B-. 

In the perturbative RG approach the normal-extraordinary equivalence 
manifests itself as follows: If hi o > 0, then its renormalized analog, hi, is 
driven to +00 under the RG flow. The variable c, on the other hand, tends to 
±00 or stays at zero, depending on its initial value. It turns out that the free 
propagator and the zero-loop profile m(z) for all these cases with h\ = 00 are 
given by the same expressions and identical with their analogs for h\ = and 
c = —00. Hence the respective perturbation series agree to arbitrary order. 
This shows that the limiting probability distributions are the same; in this 
sense there is just one fixed point for both transitions. 

The asymptotic behavior of the order-parameter profile can be written in 
the scaling form Jl8|), but with different scaling functions -P±(C) for T —* ± . 
These functions have been computed by meana_of RG-improved perturbation 
theory to one-loop order in 4 — e dimensionsEil The results, extrapolated to 
d = 3, are in reasonable agreement with Monte Carlo resultsEj The short- 



9 The reader should be cautioned that the physically reasonable distinction between the 
extraordinary and normal transitions is frequently not made in the literature: It is a com- 
mon, but unfortunate, practice to use the name "extraordinary" for both the normal and 
extraordinary transitions as well as for the associated universality class. 
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distance behavior of the scaling functions, 



P±(0 « C P/ " c± + ai,± C 1/v + aa,± C 2/ " + &± C d + ■ • • , (22) 

is in conformity with the r dependence of mi given in Eq. (p0|). The first 
three contributions correspond to the terms regular in r listed in Eq. ( pl[ ) ; the 
fourth one contains the |r| 2 ~ a singularity. This latter term may be understood 
as the contribution from the component T zz of the stress energy tensor to the 
boundary operator expansiom'E3 (BOE)0 

^..^CiWOjta) (23) 

of <p for x — (xu,z) — > a; 8 = (xm,0)@ Here (4>(x)) c oc z~^l v is the critical 
profile. The sum runs over a complete set of boundary operators with scal- 
ing exponents A\. The functions C\(z) vary cx z Aa at criticality. Consistency 
with Eq. (E2) requires that the boundary operator in (^3|) with the smallest 
scaling dimension, except the one operator 1, has A\ = d. This is T zz . 

To my knowledge, there is still no experimental system that has been 
clearly demonstrated to have supercritical surface enhancement. Thus it is 
fortunate that extraordinary critical behavior can be seen at normal transi- 
tions. A much studied phenomenon, which occurs at such transitions, is the 
critical adsorption of fluidsr^Tt^ This occurs when, for example, a binary 
liquid mixture is brought to its bulk critical point in the presence of an ex- 
ternal wall (e.g., container wall) or other distinct physical interface. For such 
a mixture, the order parameter is a composition variable. The wall usually 
favors one of the components, a property which translates into the presence 
of a surface ordering field /ii.o- Hence one gets back to our previous choice 
of Hamiltonian for the normal transition. One important signature of critical 
adsorption can be read off from Eq. (|22): the composition varies as z~@l v for 
z <C £■ A second one is that the excess order parameter m s (corresponding 
to the "total amount of adsorbed order") diverges as r~( !y ~' 3 ). The analyti- 
cal and Monte Carlo results given in Refs. @ and @ have proven useful 



for analyzing experiments on critical adsorption. For example, values for the 
universal amplitude c+ in Eq. (|2^) and for certain universal ratios involving 
integrals of P± have been extracted frjom_.experiments. These are in fair agree- 
ment with the theoretical prediction£3'Eil considering the still large error bars 
of both types of estimates. 

h The operators here should be read as renormalized ones. In the sequel we will use the label 
'ren' to distinguish renormalized from bare operators when necessary. 
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4 Field theory, boundary conditions, and short-distance singulari- 
ties 



Since a detailed exposition of the field-theory approach to surface critical be- 
havior may be found Ref. ^, I will restrict myself here to giving a brief summary 
of some essential points for readers with a background in field theory. As we 
have seen, 'operators' such as <j>{x) generally have different scaling dimensions 
for points x in the interior, int Ai = M. \ B, and on the boundary, B. For this 
reason we introduce separate bulk and boundary source terms, defining 



Z[J, Ji] = e g[J ' Jl] cx / £><•/> exp -H+ [ dV J<j)+ [ dAJ x 



(24) 



The functional Q generates the connected correlation functions G^ N ' M ^ = 
((/>... 4>4>b ■■■ 4>b) con of N fields 4>{xi) with x { e mtM and M fields 4>b = 
4>(x®) on the boundary. 

Exploiting the invariance of Z under changes (f> — > <p + 5(f> in a standard 
fashion gives the 'equations of motion' 

J(x) = -A<t> + U'(ct>) , xe'mtM, (25) 

and 

J 1 {x B ) = -d n 4> + U l 1 {4> B ) = -d^ + cafo-hw, x B eB, (26) 

where d n is the derivative along the inward normal. The resulting boundary 
condition for the zero-loop profile m'°l reads [d z — Co] m^(z = 0) = — /ii,o! the 
one for the free propagator G = [— A + U"(m^)]^ 1 (x,x') is of the Robin type, 

d n G(x B ,x') =U'{{mW)G{x B lX ') =c G{x B ,x') . (27) 

It is well known that G = Gb + G s , where Gb{\x — x'\) is its bulk analog 
while G s {x : x'j is an image term. The former has an integrable singularity at 
x = x' , giving rise to the familiar ultraviolet (uv) bulk singularities in Feyn- 
man graphs. The latter also becomes uv singular, but only for x = x' 6 B. 
This causes additional uv singularities in Feynman graphs. Owing to the local 
form of the primitive divergencies, they can be absorbed by local boundary 
counterterms . The upshot is that the required counterterms can be written 
as a sum J M CTfc(a;) + J B CTi(x) of bulk and boundary contributions. Pro- 
vided one can convince oneself (e.g., by power counting) that they all have 
the form of the interaction terms included in the action, the theory is renor- 
malizable. In particular, our model defined by Eqs. (|s|)-(|T^) turns out to be 
(super-)renormalizable for d = 4 (d < 4). One important difference with the 
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usual rcnormalization of infinite-space models should be stressed._bfjw.ever: in 



general, one-particle reducible renormalization parts may occurut 

The counterterms needed to renormalize G^ N,M ^ correspond to reparame- 
trizations of the form 

4>=Zy 2 4> Tm , to =t c + [i 2 Z t t , u = ^Z u u, (28) 

<t> B = [Z^] 1 / 2 (0 B ) ren , c = c sp + f iZ c c, h h0 =fx d / 2 (Z !t) Z 1 )-V 2 hi, (29) 

where \i is the momentum scale. The Z factors are meromorphic in e, provided 
dimensional regularization is used. In a theory regulated by a momentum cut- 
off A, one has c sp ~ A, analogous to the familiar behavior r c ~ A 2 . For the 
dimcnsionally regulated theory in fixed dimensions al < 4, a nonperturbative 
shift of the form c sp = u^ e C(e) occurs i3 comparable to Symanzik's mass- 
shifHT c = U o /e r(e). 

RG equations for the renormalized functions G^n^ follow in a standard 
fashion and can be utilized to derive their scaling forms and the scaling relations 
for the exponents of the special transition. The scaling functions involve the 
scaling variables c|t|~* and /ii|r|" A i . In addition to the bulk correlation 

length £, one has two surface-related lengths £ c ~ c~ u ^ tl? and £i ~ h 1 v ^ Al 
that become arbitrarily large as the multicritical point Sp is approached. 

For c > and h\ = 0, a crossover to a behavior characteristic of the 
ordinary transition occurs. The information about the latter is contained in 
the above-mentioned scaling functions. However, it cannot be deduced just from 
the RG equations of the Gi^l ■ To extract it, the limiting behavior of their 
scaling functions for c|r|~* — * oo (£ c <C £) must be known. Setting c to its 
value c* rd = oo at V* ld does not really help because the bare and renormalized 
theories satisfy the Dirichlet boundary condition (f>g = for cq = oo and 
c = oo , respectively, so that all functions G^J^f with M > vanish in this 
case. Hence, to determine the critical behavior of these functions for c > 0, 
one must move away from V* ld , considering large but finite values of c. An 
alternative strategy is to infer the thermal singularities of local properties on 
B from those of their analogs for small z > 0, with c = oo. This involves 
knowledge about the corresponding short-distance singularities. 

Fortunately, a ccuxvenient way of obtaining the required information has 
been found long agooO Using the boundary condition (p6]) [or (|27|), in pertur- 
bation theory], the operators <f>s i n G^ N ' M ' may be replaced by c ~ 1 dn4>- It 
follows that the desired critical behavior is given by the functions G^' M ' 1 = 
(</>... <j)d n (j) . . . d n (f>) con with c = oo. 
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The Dirichlct boundary condition of the regularized theory for cq = oo 
holds even if hi o =/= 0, provided /ii,o/co — ► as Co — > oo. In order to retain 
a nonvanishing relevant surface field g\ ~ /ii.Oj the limit cq — > oo should be 
taken with hi tOC = % hi t o fixed. Noting that the term oc in TL can be 
rewritten as —hi >00 J g d n (f>, one sees that the relevant surface operator to which 
<7i ~ fti,oo couples is <3 n 0. 

Renormalization of the cq = oo functions G^' M ^ can be achieved by means 
of the reparametrizations^ 

9„0 = [^Zl,oo] 1/a (a„0) ren , /ll.oo = M (d " 2)/2 (^^l,oo)- 1/2 /l r i:Sc (30) 

and those given in Eq. J2jf) . The implied RG equations of the renormalized 
functions G^i-cr? can be exploited ta-ifte usual manner to analyze the critical 
behavior at the ordinary transitionptrcl In particular, the scaling relations for 
the critical exponents, and scaling forms such as ( [l6| ) follow, with h 1 ^ playing 
the role of gi . 

We close with a brief discussion of short-distance singularities (cf. pp. 
190-202 of Ref. |^). Consider first the case of the special transition. We are 
interested in the asymptotic z dependence of the functions Gren at criticality, 
hi = c = t = 0, as one of the N off-surface points, x = (acn,z), approaches 
x B = (sen, 0). The answer is provided by the BOE 

^ cn (x) » C{z) {0 B Y cn (x B ) with C{z) oc z -W-fin/ v . (31) 

z^O 



The form of C{z) follows from its RG equation, implied by those of G^l . 
The exponent of z can be rewritten as (/3i — (3)/v = r)± — r/, using scaling laws. 
The analog of this BOE for the functions G^i-an at = r = readJ 

P cn (x) « 0^(8^)^ (x B ) with Goo (z) cx z^-W- . (32) 

z^0 



Obviously, this BOE applies equally well to the Gren at c = oo and hi = 
t = 0. That the leading contribution now arises from d n <f> rather than from <^g 
reflects the Dirichlet boundary condition of the renormalized theory at c = oo. 

As one moves away from the fixed points "P* and "P* rd , permitting some 
relevant fields to be nonzero, the coefficient functions appearing in the operator 

'For \c\ < oo one can easily see from the boundary condition that a renormalized operator 
(dn<f>) rca can be defined by d n 4> = Z c [Z t f,Zi] 1 ^ 2 (d n 4>) rca , in the dimensionally regularized 
theory. However, for c = oo this is no longer true. 
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algebra, in general, also depend on these scaling fields. Further, additional 
contributions may occur. Let us specifically consider what happens to the 
short-distance singularities of the functions Got^'ren and the BOE ([32|) when 
a small surface held, is turned on. This problem was investigated and 
solved by Symanzhc3 in 1981. His results do not seem to have been recognized 
widely, especially not in the solid state community, even though they were 
reformulated in the language of the BO.E and generalized to manifolds with 
curved boundaries in subsequent workliU 

Note, first, that the analog of Eq. ( |26| ) for the cq = oo functions Go^rer? 
reads 

Ji,oo(x B ) = (j)® - hi t00 , (33) 

where Ji.oo is the source associated with d n cj). In the absence of Ji,oo, this 
reduces to the boundary condition 0g = hx t00 for the regularized bare theory. 
(Position-dependent boundary values 4>b can be imposed by retaining a nonzero 
</i,oo, as was actually done by Symanzik.) The naive analog of this boundary 
condition for the renormalized theory, limaj^^e </> lcn = h^^, does not hold 
because of short-distance singularities. Instead one has (setting = 1) 

Ctiz)- 1 <f Ten (x) — h\% , (34) 

2^0 



where the choice Ct(z) = J B dAGx^ln(x,0) with h 1 ^^, = can be madej-Zc 
pp. 16-17 of Ref. ^5|). The behavior of this function at V* vd is well knowitrer 



and follows directly from its RG equation, explicitly given in Ref. 25. One 
finds 

C t (z) oc z 1 -^ . (35) 

The exponent 1 — r]± d = (A° rd — (3)/v reflects the different asymptotic scale 
dependencies ~ £- A 7 d / l/ and ~ 1&I V of K^oo ano - ^ ren m the infrared limit 

Thus the BOE become£3 

<t>™(x) « C % {z)h\% +C 00 (z)(d n cf>y™(x B ) + ... , (36) 

2^0 



where the ellipsis stands for less singular contributions (arising from other 
boundary operators as well as from the omitted additional dependence of co- 
efficient functions such as C\ and C^ on the scaling variable W-y^ z Al .) 

The existence of the short-distance singularity (|3^) is a central issue of a 
recent Letter|23 in which it is verified by means of Monte Carlo simulations for 
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the d = 3 Ising_£ase. In subsequent work by the same authors the d = 2 Ising 
case is studiedE3 

The information acquired about short-distance singularities can be utilized 
to predict the position dependence of the order-parameter profile m(z) = {<j>) 
in various asymptotic regimes with z <C £. At hulk criticality, r = 0, the 
scaling form of the profile, m = m c , reduces tcEltH 

m(z) « z~^ v F(zh v '^ \chf /AT ) . (37) 

The limits F(oo, z-i) and F{z\, 0) of the scaling function F[z\, Z2) must exist 

and be nonzero. Further, we know from the BOE ( |3l| ) that F(zi,0) ~ zf 1 ^ 
as zi — > 0.^] Consistency with the scaling form near V* ld requires that 

F( Zl = zh^,z 2 = ch^) ~ W(zg^) (38) 



wit 



li 



, 9l =^c-, tf =^L^ = m_J*L, (39) 



where we know from Eqs. (35) and (|36|) that W(z^) ~ z\ n± as 2:3 — > 0. In the 



scaling regime where Eq. (37) holds (small c and hi), both lengths £ c ~ c 



-u/9 

and £1 ~ h x are large. Eq. ( fj8| ) involves a length £ or d ~ ff^"^ 1 i which 

can be large or small. The resulting power laws describing the z dependence 
of m c {z) in the various regimes are 

\z-W-f*F>l" forz«£ c , fc, £ ord , 
m c (z) - ^ zi-^ for £ c < z < £ ord , (40) 

I for i c « £ ord « z . 

These findings are in conformity with unpublished one-loop results for 
m c (z), obtained some time ago by Ciach and myseliS Details can be found in 



a 



recent paper by Ciach and RitschelE3 



5 Concluding remarks 

As I have tried to illustrate in this talk by discussing prototypical examples 
of such phenomena, the application of field-theoretical RG methods to the 
study of boundary critical phenomena is well established. It has led to deep 

^Notc that in Eq. (3.227) of Ref. ^jthe exponent of this power contains an incorrect minus 
sign; at other places of this reference the correct value is given. 
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insights, has provided a reliable theoretical framework, and yielded predictions 
for experimentally accessible quantities. Research over the past two decades 
has revealed a surprising wealth of such phenomena. 

It is particularly encouraging to see that the number of experiments de- 
voted to careful investigations of boundary critical phenomena has been con- 
stantly increasing for some years. Since more precise experimental data should 
be available soon, an obvious challenge for theorists is to improve the still rather 
moderate accuracy of the field-theoretical predictions. To this end, it would 
be very desirable to extend the presently available two-loop— results for the e 
expansions and the massive RG approach in fixed dimensional d < 4 to higher 
orders, so that more precise field-theoretical estimates for critical exponents 
and amplitude ratios can be gained via Pade-Borel resummation techniques. 

While achieving greater accuracy is certainly one important goal for fu- 
ture research, it should also be emphasized that only relatively few models for 
boundary critical phenomena have been thoroughly investigated so far. Thus 
interesting new features may still have to be discovered. 

To illustrate this point, let me finish by reporting a result established 
recentlyEaO As we have seen above, systems exhibiting boundary critical be- 
havior at bulk critical points also can be divided into universality classes. 
These surface universality classes depend on the bulk universality class and 
additional relevant surface properties. Examples of the latter we have come to 
know so far are: the short range of the change of interactions caused by the 
surface; their short-range nature; whether the surface enhancement is subcrit- 
ical, critical, or supercritical; and whether symmetry-breaking surface terms 
(with h\fi ^ 0) are present or not. 

In the case of antiferromagnets in a magnetic field (sufficiently weak so 
that the bulk transition remains continuous) and binary alloys with non-ideal 
stoichiometry and a continuous order-disorder bulk transition, it turns out that 
the surface universality class in general also depends on the orientation of the 
surface plane with respect to the crystal axes. There exist both orientations 
that preserve, as well asJJiose that break, the symmetry with respect to the 
two types of sublatticesu Depending on whether the orientation is symme- 
try preserving or breaking, ordinary and normal transitions (corresponding to 
the ordinary and extraordinary universality classes) are expected to occur, re- 
spectivelxlrjrovided the enhancement is subcritical, which seems reasonable to 
assume)E3lI3 To appreciate this result one should note that, on the level of a lat- 
tice description, no ordering bulk and surface fields are present, because these 
would correspond to staggered magnetic fields, in antiferromagnetic language. 
It should be possible to check these predictions by means of experimental stud- 
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ics of the continuous A2-B2 bulk transitions in binary alloys such as FeCo[] or 
FeAl. 
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Note added in proof 

We mentioned above the result taken from Ref. El that the disorder-order tran- 
sition of FeAl from the B2 to the DO3 phase should belong to the universality 
class of the 0(2) |</>| 4 model. Recently R. Leidl (U. Essen, unpublished) has 
reanalyzed the problem within Landau theory. He found that the Landau 
theory allows quadratic anisotropics. These are relevant and should drive the 
system to the the n = 1 (Ising) fixed point. Thus the transition should belong 
to the universality class of the one-component |<?!>| 4 model. Upon contacting 
Prof. Binder's group we learned that similar conclusions have been reached in 
the paper by W. Helbing, B. Diinweg, K. Binder and D. P. Landau, Z. Physik 
B 80, 401 (1990), which corrects an earlier one by B. Diinweg and K. Binder, 
Phys. Rev. B 36, 6935 (1987) whose different conclusions are described in Ref. 
48. We are grateful to Prof. Binder and B. Diinweg for kindly informing us 
about this matter. 
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